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EQUATION OF THE FREE-FALLING NAPPE 


Fred W. Biaisdell,! M. ASCE 


SYNOPSIS 


A general equation for the form of the nappe is developed in this 
paper. The equation is not valid close to the crest, but applies to that 
portion of the nappe that is free-falling—where pressures within the 
nappe are atmospheric. The constants in the equation have been evalu- 
ated for the vertical, sharp-crested weir having approach channel 
depths ranging from deep to zero (the free overfall); that is, over the 
entire range of subcritical approach velocities. The evaluation is made 
using nappe coordinates obtained by others. The equation is checked 
by comparing its predictions with a number of published profiles. The 
comparison is shown to be excellent for the lower approach velocities 
and good at the higher approach velocities. 


INTRODUCTION 


An adequate number of carefully performed experiments already 
have been made and published which give the form of the nappe pro- 
duced by a vertical sharp-crested weir. It would seem that any further 
writings on the subject would be superfluous. However, no general 
equation for the form of the nappe that covers velocities of approach 
other than zero has been derived from the experimental data. 

The need for a general equation for the form of the nappe became 
apparent when the writer and his colleagues were assigned the task of 
developing a design for a stilling basin for use with straight drop spill- 
ways. This is because it is desirable to use general formulas for the 
determination of the component dimensions when developing general- 
ized designs. The straight drop spillway is ordinarily a vertical weir 
over which water falls onto an apron on which energy is dissipated. 
The point at which the nappe strikes the apron is important since this 
in some measure determines the location of the blocks and sills used 
to assist the energy dissipation. It is also a factor in the determination 
of the basin length. 


1. Project Supervisor, U. S. Department of Agriculture, ARS, St. 
Anthony Falls Hydraulic Laboratory, Minneapolis, Minnesota. 
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Scope 


The study reported here has resulted in the development of a general 
formula for the trajectory of the lower and upper nappes of a vertical 
sharp-crested weir with velocity of approach (or approach channel 
depth) as a variable. The equation applies only to that portion of the 
nappe which is free-falling, that is, where the pressure at every point 
within the nappe is atmospheric. The equation is valid for approach 
velocities ranging from zero to the critical velocity. 

It is well known that convergence of the streamlines causes the 
pressure within the nappe to be greater than atmospheric for a short 
distance downstream of the crest. The formula presented here repre- 
sents a parabolic curve. This curve form is produced when the only 
force acting on the nappe is that of gravity. Therefore, the formula is 
not valid in the vicinity of the weir crest where pressures within the 
nappe are above atmospheric. The lower limit of the applicability of 
the equation is at x/H=0.5. Here x is the abscissa measured from 
the upstream face of the weir and H is the specific head. 


Works Utilized 


Data obtained from a number of different sources were analyzed by 
the writer. 

Probably the best known nappe coordinates that are used here are 
those which William P. Creager (1),2 Hon. M. ASCE, extrapolated from 
Bazin’s experiments. Bazin’s data extended only about to x/H = 0.12 
for the upper nappe and to about x/H = 0.65 for the lower nappe. This 
is probably within the range of positive pressures within the nappe. 
Mr. Creager extended the nappe by considering the horizontal velocity 
to be constant and the vertical velocity to increase under the influence 
of gravity. His lower nappe formula for the vertical weir with zero 
approach may be written 


2 
= 0.190 + 0.046(5) - 0.366(=) (1) 


where y is the ordinate measured from the weir crest. Positive val- 
ues of y are above the crest level, negative values below. 

The Bureau of Reclamation (2) of the U. S. Department of the Interior 
conducted a most thorough and detailed study of the nappe form. Their 
data were analyzed by plotting the nappe coordinates to a large scale, 
drawing a smooth curve through the data, reading coordinates from the 
curve and dividing by H, combining the coordinates according to h,/H, 


plotting the coordinates for each value of h,/H, drawing an average 
curve through this data, reading coordinates from this average curve, 


2. Numerals in parentheses, thus: (1), refer to corresponding items 
in the Bibliography (see Appendix). 
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and tabulating the readings. Here h, is the velocity head in the ap- 


proach channel. The data given in Tables 12 and 18 of Reference (2) 
were used by the writer. 

Hinds, Creager, and Justin (3) present a table of nappe coordinates 
for various approach velocities. Their table was adapted from a thesis 
by Jacob E. Warnock (4), M. ASCE. Mr. Warnock’s thesis was based 
on the Bureau of Reclamation data mentioned previously. The close 
agreement between the two sets of data that later will become apparent 
is therefore not surprising. 

Calvin V. Davis (5), M. ASCE, also presents a table of nappe co- 
ordinates for vertical sharp-crested weirs. The data on which this 
table is based were obtained by the Bureau of Reclamation, M. Bazin, 
and Professor Scimemi. In a letter to the writer, Mr. E. W. Lane, 

M. ASCE, who prepared the table, states that the data he used were not 
the Bureau of Reclamation data mentioned previously. It was data ob- 
tained at an earlier date. It may, therefore, be considered as distinct 
from the other Bureau of Reclamation data. 

Ettore Scimemi (6), M. ASCE, has published experimental data ob- 
tained on a vertical sharp-crested weir. The particular value of 
Professor Scimemi’s presentation is that his nappe coordinates cover 
a much wider range than does any other data that the writer has been 
able to discover. For the lower nappe, Professor Scimemi gives an 
equation which can be written 


2 
= 0.182 - 0.054(=) - 0.1416(=) for = > 0-50 (2a) 


while for the upper nappe 


2 
% = -0.243 + 0.694(2) - 0.496(2) for >1.h0 (av) 


H H 


A. T. Ippen (7), M. ASCE, shows the form of the nappe in his Fig. 11, 
p. 529 (this figure was originally presented by Hunter Rouse (8), M. 
ASCE) and gives a chart of nappe coordinates in his Figs. 12, pp. 530- 
531. Fig. 11 gives the nappe form for four approach channel depths. 
The coordinates of these nappes were scaled from the figure. In addi- 
tion, coordinates of the nappe were read from Figs. 12 for a number 
of approach channel depths which were subsequently converted to ap- 
proach velocity heads. The data obtained from each of these figures 
were used by the writer. 

The final group of data used by the writer was unpublished informa- 
tion obtained by Donald A. Parsons, a former colleague of the writer. 
Mr. Parsons’ data apply to the free overfall only. 


Method of Analysis 


In common with Mr. Creager, Dr. Ippen, and others, the writer 
makes the reasonable assumptions that the horizontal velocity is 
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constant and that the only force acting on the nappe is gravity g. A 
particle of water in the lower nappe will, in the time t, travela 
horizontal distance x from the face of the weir of 


x =v, t cos @ (3) 


In the same time the particle will travel a vertical distance y of 


+v 


In Eqs. 3 and4 v a is the velocity, 6 is the angle of inclination of the 


velocity with the horizontal at the point where x = 0, and C' is the 
value of y at x= 0. Solving Eq. 3 for t, substituting in Eq. 4, and 
simplifying 


y = C' +x tan © + g x (5) 


cos’: 
° 


Eq. 5 gives the form of the lower nappe. 

Since the horizontal velocity is constant the vertical thickness of 
the nappe T must also be constant. The equation of the upper nappe 
is then 

2 
y = C! +x tan @ = +T 
2 ~ (6) 


It is advantageous when developing a general formula to use dimen- 
sionless units. The basic unit of measurement will be the specific 
head H. This is the total head over the weir—the observed head plus 
the velocity head. Dividing both sides of Eqs. 5 and 6 by H and sub- 
stituting C for .C'/H, and S for gH/v,” cos? @ the equation for the 


lower nappe is 
2 
Ze (x 
tan +5 (5) (7) 


and for the upper nappe 


2 
S (x 
C +5 ten + (5) + 


Eqs. 7 and 8 are the basic equations used in this analysis. 


3. It is not necessary to specify any particular part of the nappe because 
the development is completely general. The lower nappe is specified 
here because its form is desired initially. 
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The method of analyzing the data was suggested by Donald A. Parsons. 
It is so simple and such a useful tool that it is surprising that apparently 
no previous use has been made of it, at least so far as the writer can 
deduct from the available publications. 

Differentiating Eq. 7 and rearranging 


d(y/H x 
(= 9 
dix = tan @ (9) 
An identical result is obtained on differentiating Eq. 8. This indicates : 


that coordinates of both the upper and lower nappes may be used to 
evaluate tan @ and S. Steps in the analysis included: (1) reducing ob- 
served values of x and y to x/H and y/H; (2) computing A(x/H), 


A(y/H), and ace. (3) plotting se vs. x/H as in Fig. 1; > 


fitting a curve (a straight line) to the plotted points; and (4) determining 
S from the slope of the curve and tan 6 from the intercept of the curve 
at x/H=0. The constant C was computed for each plotted point by 
substituting tan@ and S in Eq. 7. 
The scale of the plots made it difficult to read tan@ and S as 
precisely as was desired. Therefore, C was plotted against x/H and 
(x/H)? to check tan@, S, and C. A typical plot for the data appear- 
ing in Fig. 1. is shown in Fig. 2. The slope in Fig. 2a indicates that 
tan @ was read incorrectly and that its true value is 0.085 instead of 
0.08 as was read originally. Also from Fig. 2a the correct value of C 
is seen to be 0.139. If the data in Fig. 2a had curved systematically it , 
would have indicated anerror in S. The plot of Fig. 2b has been pre- i 
pared to indicate that S is satisfactory after correcting tan 6 . The 
waviness of the data in Figs. 1 and 2 is an indication that the Bureau of 
Reclamation table, in spite of the care lavished on it, is still somewhat 
rough. Mr. Lane smoothed his table by computing second differences 
and his results do not exhibit the scatter shown in Figs. 1 and 2. 
The nappe thickness T was computed by substracting the upper 
nappe coordinates from the lower nappe coordinates and averaging the 
results. The variation from the average was not great until the nappe 
surface appeared to break up and accurate surface readings became ' 
difficult. The average nappe thickness between x/H values of 0.5 and 
1.8 is the figure recorded for the Bureau of Reclamation and the Hinds, 
Creager, Justin data. For the other data employed, the nappe thickness 
-was computed for values of x/H > 0.5. 
The equations given in this paper are valid only at values of 
x/H > 0.5. This is a conservative general statement. The plotted data 
indicate that the equation may be valid at lower values of x/H for the 
lower nappe than for the upper nappe. 
The reason why the equation is not valid between x/H values of 0 
and about 0.5 is apparent after considering the difference between the 
original assumption of atmospheric pressure within the nappe and the 
actual condition existing near the crest. Professor Scimemi’s (6) “7 
experimental results show that pressures within the nappe are greater 
than atmospheric to 0.5x/H but that at about this point they become 
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negligible and that they disappear at higher values of x/H. Therefore, 
it is not to be expected that the equations developed here will be valid 
close to the weir crest. 

A physical illustration of the meaning of the constants in Eq. 7 will 
perhaps be helpful. The constants C ard tan @ are illustrated in 
Fig. 3a. There it can be seen that C is the intercept of y/H at the 
plane of the weir crest. The value of C would be zero if the nappe 
were free-falling at the weir crest. The nappe is under the influence 
of converging streamlines at the crest and is not free-falling. Thus C 
has a value other than zero except for the case of the free overfall 
where C approaches zero. The constant C is found, as was men- 
tioned previously, by substituting in Eq. 7 after tan @ and S had been 
evaluated. Also the controlling streamline at the plane of the weir 
crest is not horizontal. For the vertical weir, it is inclined slightly 
upward while for the free overfall it is inclined downward. The in- 
clination is illustrated in Fig. 3a and its magnitude is found by pro- 
jecting the slope of the free-falling nappes, shown in Fig. 3b, back to 
the plane of the weir crest. The slope, tan 6, can also be found by 
substituting in Eq. 9 after S has been evaluated. 

The vector diagram of Fig. 3c explains S. The velocity vectors at 
points 1 and 2 in Fig. 3a are transferred to Fig. 3c. The change in 
velocity direction between points 1 and 2 is a result of the influence of 
gravity acting on the vertical component of the velocity, since an origi- 
nal assumption was that the horizontal component of the velocity was 
constant. The vertical acceleration, represented by S, was actually 
determined from the slope of the curve of Fig. 3b. It could have been 


determined also by computing the second differential of Eqs. 7 or 8 
which is 


2 


a H 
(10) 


d (x/H) 
Results of Analysis 


Form of the Nappe for Zero Velocity of Approach 


More data are available for the vertical weir having a zero velocity 
of approach than for other approach velocities. It is therefore possible 
to provide a more reliable determination of the constants in Eq. 7 for 
this condition than for other velocities of approach. The original data 
for each of the lower nappe profiles were first plotted so they might 
be compared. These are shown in Fig. 4. It will be noted that the pro- 
files compare well except for the Davis (Reference (5), p. 338, Table 1) 
and the Creager (Reference (1), Fig. 25) profiles. The Davis curve 
falls below the other curves. Since the Bureau of Reclamation data on 
which this curve is largely based were superseded by data obtained 
later by the same Bureau, the Davis coordinates were not used in 
computing values of the constants. The Creager data fall above the 
other curves. Since these coordinates were projected by Creager 
from data by Bazin that extended only to x/H ~ 0.65, they too were 
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not averaged. In this connection, Mr. Creager (1), p. 110, notes, “The 
paths are only approximate, as they are extended from experimental 
points relatively close to the top of the dam. . .” 

The values of C, tan 9, S/2, and T determined from the various 
data are given in Table 1. The averages for the Scimemi; Hinds, 
Creager, and Justin; Bureau of Reclamation; and Ippen data ave also 
given. Also listed are values of the constants which the author proposes 
as being even values that are close to the average. The curve in Fig. 4 
is plotted using these proposed values. The excellent agreement with 
the observed data is obvious. 


Effect of Approach Velocity 


Increasing the velocity of approach causes reductions in C, in tan 6, 
and in the effect of gravity. The effect is shown in Fig. 5 where C, 
tan 6, S/2, and T are plotted against the relative velocity head 
h,/H. Only the Bureau of Reclamation; Hinds, Creager, and Justin; 


and Ippen give nappe coordinates for different approach velocities. The 
trend of the data was established by the plotted points. The initial point 
on each curve was established by the proposed value given in Table 1. 
The equations for the four constants in Eqs. 7 and 8 are 


h 


= 0.150 = (1) 


h 
tan @ = 0.411 - 1.603(-=) 


h 2 h 
- / 1.568 (+7) 0.892(-) + 0.127 


S hy 
= -0.4425 + 0+25(- 


(14) 


where m= (h,/H) - 0.208. Substituting these values in Eqs. 7 and 8 
gives the general equations for the lower and upper nappes for the 
vertical weir. The equations are, for the lower nappe, 


\ 
and : 
482-7 


ho 
= 0.150 0.45(+F) + | 1.603(-7) 


= 0.892(--) + 0.127 | 


bh, x 2 
+ | -0.425 + | (=) 


and, for the upper nappe, 


pa hy 
= 0-150 - 0.45(--) + | 0.411 - 1.603(-5 ) 


x 
- 1.568(—-) - 0.892(—-) + 0.127 


2 10m 


h 2 
+ 1:25 + | =) > 0.57 


(16) 


These equations are valid only for x/H>0.50 and for values of h,/H 
between 0 and 1/3—the subcritical velocity range of flow. 


Comparison with Published Nappes 

The proof of Eqs. 15 and 16 is in the comparison of their predictions 
with the published nappe profiles. The comparison for h,/H =0 is 
shown in Fig. 4 to be excellent. Comparisons for h./H = 0.10, 0.20, 
and 0.33 are shown in Fig. 6. The comparisons in general are seen to 
be very good. Only at h,/H = 0.319 was the comparison poor (not 


shown). There is need for more data between relative velocity heads 
of 0.2 and 0.33 to better define the equation between these limits; the 
equation admittedly needs additional verification within this range. 


Nappe Form for 0<x/H<0.5 


It has been noted previously that Eqs. 7 and 8 or Eqs. 15 and 16 are 
not valid when x/H is less than about 0.5. If the form of the nappe is 
needed below x/H = 0.5 it is necessary to obtain it from some other 
source. It is suggested that the coordinates be obtained from one of 
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the sources that have shown good agreement with the equations in Figs. 
4 and 6. Some slight adjustment may be required to insure that the two 
curves join perfectly, but no adjustments that would cause any concern 
are anticipated. 


SUMMARY 


Eqs. 15 and 16 give the form of the nappe which falls freely from a 
vertical, sharp-crested weir. The nappe is free-falling beyond about 
0.5H from the weir crest. The equations are valid for approach 
velocities lying within the entire subcritical velocity range; they are 
valid for approach channel depths ranging from deep to zero below the 
weir crest. 

The equations give excellent agreement with experimental data for 
the lower approach velocities and good agreement for all but one of the 
higher approach velocities. 


COMMENT 


The fact that the equations are not valid in the vicinity of the crest 
may be considered a serious deficiency by some persons. It should be 
remembered that they were not developed to give the form of a dam 
face, although Eq. 15 can be used for this purpose. The equations were 
developed to give the point at which the free-falling nappe hits the apron 
in order to provide one of the dimensions necessary for the determina- 
tion of the apron length. 
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TABLE I 


COMPARISON OF NAPPE EQUAT IONS 


= 0 


Reference (1), p. 108 


Reference (1), p. 110, 
Table XVIII 


Reference (2), Tables 12 & 18 
Reference (3), p. 363 
Reference (5), p. 338 
Reference (6), H = 4. om 
Reference (6), H = 8.8 om 
Reference (6), H = 13.2 om 
Reference (7), p.- 529, Fig. 1l 


Reference (7), pp. 530-531, 
Fig. 12 


Averages of References 


(2), (3), (6), (7) 
Author 


0.173 


0.190 


0.139 
0.153 
0.02 
0.130 
0.16 
0-16 
0.198 
0.190 


0.157 


0.610 
0.588 


0.538 


Source tan @ 8/2 
MM +0.057 -0.3675 
40.085 -0.1,0 
+0.060 -0.430 -- 
+0.090 0.578 
+0.050 -0.422 0,593 
+0.050 0.568 
40.050 -O.410 0.512 
0.030 =0.405 0.573 
-0.035 -0.418 0.560 
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(a) Test of tan © 
e tan 6 = 0.080 x tan 6 = 0.085 


(b) Test of S 
$+ 0.86 


Fig. 2— Tests of S and tan © 


| 
e Lower nappe 
| 
4 50 2 3 
x/H 
a Fig. |— Typical Plot for the Determination of S and tan 6 
0.17 
x/H 
| 
0.17 
q 
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(a) Section Through Nappe 


44/4) 
ax, 7 


vet = A(x/H) 


1 /H) 


ton" 


diyg/H) _ 
d(x,/H) d(x,/H) 


A(x/H) S A(x/H) 


(b) Graphical Determination of tan ©@ and S (c) Vector Diagram 


Fig. 3— Illustration of C, ton 6, and S 


= 
7 
: 
A tan © x/H 
= 
(y/H) 
tan x/H wit 
: 
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Fig. 25 

Table i2 

3], Page 362, Table | 
5], Page 338, Table | 
[6], Table 5, H= 4.40 cm 
[7], Page 529, Fig. II 
[7], Page 531, Fig. 
D. A. Parsons 


© 
° 
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v 
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x/H } 
2 3 a 5 6 7 
| 
| 
1% Fig. 4— Lower Nappe Comparison for h,/H=0.00 | 
15 
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Fig. 25 
(2], Table 12 
[3], Page 362, Table | 
[5], Page 338, Table! 


[6], Table 5, H= 4.40 cm 
[7], Poge 529, Fig. 11 
[7], Page 531, Fig. i2b 
D.A.Parsons 


ig. 5— Effect of Approach Velocity on C, tan ®, S/2, and T 


cas 
0.20 0.15 
0.16 by 0.10 
a 
0.12 0.05 
> + 
No 
0.04 0.05 
0.04 0.15 \ 
0.060 0.1 02 O38 04 0.200 0.1 0.2 0.3 ©0.4 | 
h,/H by/H 
0.22 0.68 
+ 
0.30 0.60 
0.34 ty T 0.56 va 
0.86 we ° 
02 #038 04 0.440 0.1 02. #03 0.4 
hy/H hy/H 
Reference Symbol Reference 
bd 
Vv o 
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Reference 
[i], Fig. 25 
[2], Table 12 
[3], Poge 362, Table | 
[5], Page 338, Table | 
[6], Table 5, H = 4.40 cm 
[7], Page 529, Fig. I! 
[7], Page 531, Fig. 
D. A. Parsons 


4 
| 
q 
7 
x/H 
2 
es Fig. 6 ~ Lower Nappe Comparison for h,/H = 0.10, 0.20, and 0.33 
| 
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300(EM)*, 301(SA)*, 302(SA)*, 303(SA)®, 304(CO)®, 305(SU)*, 306(ST)*, 307(SA)*, 308(PO)*, 309(SA)*, 310(SA)*, 311(SM)4, 
312(SA)*, 313(ST)*, 314(SA)*, 315(SM)*, 316(AT), 317(AT), 318(WW), 319(IR), 320(HW). 

NOVEMBER: 321(ST), 322(ST), 323(SM), 324(SM), 325(SM), 326(SM), 327(SM), 328(SM), 329(HW), 330(EM)*, 331(EM)*, 
332(EM)®, 333(EM)°, 334(EM), 335(SA), 336(SA), 337(SA), 338(SA), 339(SA), 340(SA), 341(SA), 342(CO), 343(ST), 344(ST), 
345(ST), 346(IR}, 347(IR), 348(CO), 349(ST), 350(HW), 351(HW), 352(SA), 353(SU), 354(HY), 355(PO), 356(CO), 357(HW), 
358(HY). 


DECEMBER: 359(AT), 360(SM), 361(HY), 362(HY), 363(SM), 364(HY), 365(HY), 366(HY), 367(SU)°, 368(WW)°, 369(IR), 
370(AT)©, 371(SM)°, 372(CO)°, 373(ST)°, 374(EM)°, 375(EM), 376(EM), 377(SA)°, 378(PO)°. 


VOLUME 80 (1954) 


JANUARY: 379(SM)°, 380(HY), 381(HY), 382(HY), 383(HY), 384(HY)°, 385(SM), 386(SM), 387(EM), 388(SA), 389(SU)°, 390(HY), 
391(IR)©, 392(SA), 393(SU), 394(AT), 395(SA)°, 396(EM)°, 397(ST)°. 


FEBRUARY: 398(IR)4, 399(SA), 400(CO)4, 401(SM)°, 402(AT)*, 403(AT)4, 404(IR) 4 405(PO)4, 406(AT)4, 407(su)4, 408(su)4, 
409(Ww)4, 410(AT)9, 411(SA)4, 412(PO)9, 413(HY)4. 


MARCH: 414(ww)4, 415(sU)4, 416(sm)9, 417(sM)4, 418(AT)4, 419(SA)4, 420(SA)4, 421(AT)4, 422(SA)4, 423(CP)4, 424(AT)4, 
425(SM)4, 426(IR)4, 427(ww)4. 


APRIL: 428(HY)°, 429(EM)°, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)©, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


JUNE: 444(SM)®, 445(SM)®, 446(ST)®, 447(ST)®, 448(ST)®, 449(ST)®, 450(ST)®, 451(ST)®, 452(SA)®, 453(SA)®, 454(SA)®, 
455(SA)®, 456(SM)®. 


JULY: 457(AT), 458(AT), 459(AT)°, 460(IR), 461(IR), 462(IR), 463(IR)©, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 475(SM), 476(SM), 477(SM), 
478(SM)°, 479(HY)°, 480(ST)°, 481(SA)°, 482(HY), 483(HY). 


a. Presented at the New York (N.Y.) Convention of the Society in October, 1953. 

b. Beginning with “Proceedings-Separate No. 290,” published in October, 1953, an automatic distribution of papers was in- 
augurated, as outlined in “Civil Engineering,” June, 1953, page 66. 

c. Discussion of several papers, grouped by Divisions. 

d. Presented at the Atlanta (Ga.) Convention of the Society in February, 1954. 

e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 
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